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ABSTRACT 

We have used the Parkes Multibeam system and the Sloan Digital Sky Survey (SDSS) to 
assemble a sample of 195 galaxies selected originally from their Hi signature to avoid biases 
against unevolved or low surface brightness objects. For each source 9 intrinsic properties are 
measured homogeneously, as well as inclination and an optical spectrum. The sample, which 
should be almost entirely free of either misidentification or confusion, includes a wide diver- 
sity of galaxies ranging from inchoate, low surface brightness dwarfs to giant spirals. Despite 
this diversity there are 5 clear correlations among their properties. They include a common dy- 
namical mass-to-light ratio within their optical radii, a correlation between surface-brightness 
and Luminosity and a common Hi surface-density. Such correlation should provide strong 
constrains on models of galaxy formation and evolution. 
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1 INTRODUCTION 

Searching for systematic correlations among the global properties 
of galaxies, analogous to the H-R diagram for stars, may offer the 
best hope of having a better understanding about their formation 
and evolution. The recent availability of large and systematic data- 
sets at several wavelengths makes this a good time to search. Such 
searches go back to the optical pioneers of galaxy exploration such 
as Hubble (1937), Holmberg (1965), de Vaucouleurs (1991) and 
Zwicky (1942). The trends and correlations they discovered, re- 
emerge very clearly in the Sloan Digital Sky Survey, as described 
by Blanton et al. (2003) in their analysis of 200.000 galaxies at 
redshifts of ~ 0. 1 . 

This paper deals with a much smaller sample of galaxies, but 
ones found by an entirely different technique, i.e. in a blind search 
for neutral Hydrogen gas at 21 -cm. The principal motivations for 
such a blind search are two. First of all, by definition, optically se- 
lected galaxies already contain many stars, which may not be the 
case in the younger, or less evolved galaxies, which may have the 
most to tell us about their formation and evolution. Apart from X- 
ray searches, which detect the very hot gas in the potential wells 
of giant ellipticals, most galaxy searches rely, either directly or in- 
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directly, on the presence of stars and so are biased against young 
or unevolved objects. Secondly, a blind Hi search offers a unique 
way round the strong optical surface brightness selection effects, 
which could disguise the very correlations one is looking for. By 
definition a galaxy must be separately both luminous enough and 
large enough to distinguish it above a sky background which, by 
galaxy standards, is quite bright. These two separate requirements 
squeeze optically selected galaxies into an extremely narrow range 
of surface brightnesses (SB) centred, for a given catalogue, on Ecat 
where: 



Ecat 



(1) 



where lap and 9ap are the minimum apparent luminosity, and the 
minimum apparent angular radius which the galaxy catalogue will 
accept (Disney 1999). The Full Width Half Maximum (FWHM) of 
the detectable range is typically only 3 magnitudes wide. Such a 
narrow stricture certainly impoverished the old photographic sur- 
veys (Disney 1976, Disney & Phillipps 1983, Dalcanton 1997). It 
might be naively thought that CCDs would have cured this stric- 
ture. Not so, because they detect galaxies that are both fainter and 
smaller than emulsions, their "Scats' s are (see equation (1)) scarcely 
any dimmer. On the other hand, the limitation with an Hl-selected 
sample is that it will miss the early type galaxies - which contain 
very little neutral gas. According to the SDSS survey late types 
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make up the large majority of all galaxies, although they may emit 
less than their fair share of light (Blanton et al. 2001). 

Blind Hi surveys face four main problems. To reach sources 
of a given column density of Hi {^hi in atoms cm~^) any survey 
must be at least sensitive enough to find them in the most favourable 
case, i.e. when they fill the beam. In that case (Minchin 2003, Dis- 
ney 2008) toils (per beam) > k/U^m (min)]^ where k is indepen- 
dent of dish-diameter D (because a larger dish projects the same 
system noise onto a smaller area of sky). Since there is a loose 
correspondence between Nhi and SB (e.g. Swaters et al. 2003). 
Mathematically: 

iV„ - 10^°-^ f ^) 10[°-^(^^-^-'l atoms cm-^ (2) 

(Disney & Banks 1997) where (^xf^) units and 7*21 

is the average SB in B mag arcsec"^, taken over the same area as 
Nh7. Since the effective SB (f-ieff ) of an exponential disc is 2.1 
magnitudes brighter than /i2i (Salpeter & Hoffman 1996) column 
density sensitivity at the Nhi < lO^^'^cm"^ level is needed to 
detect LSBGs (i.e. pe// 25.0 B mag arcsec"^) with ^ « 

0.3 (see Section 3), which in turn requires integration times per 
beam of several hundred seconds - which rendered earlier blind 
Hi surveys either too insensitive or too impractical to detect LS- 
BGs. Since this difficulty was not recognised tmtil recently, some 
earlier claims, based on small blind Hi surveys, to set rigid upper 
limits to the amount of cosmic Hi or the numbers of LSBGs in the 
Universe, must be set aside [e.g. Shostak 1977, Fisher and TuUy 
1981, Zwaanet al. 1997]. 

The coming of the Multibeam Hi detector (Staveley-Smith et 
al. 1996) was essential to carry out such blind Hi surveys. HIPASS, 
the Hi Parkes All Sky Survey (Meyer et al. 2004) used the first 
multibeam system to survey the entire Southern Sky, and the north 
up to -1-25°. More than 4000 sources were identified in the South- 
em Sky alone. The Equatorial Survey, described in this paper, was 
initially part of HIPASS, but with an accelerated search so as to ex- 
ploit the SDSS-DR2 optical data when it emerged. HIPASS sources 
are typically more than a degree apart which makes for a real chal- 
lenge in obtaining complementary optical, and other data, simply 
because, before SDSS, each source required a separate observing 
campaign. 

The strong clustering of galaxies, and of Hi galaxies in partic- 
ular, makes the identification of the source with an optical candidate 
quite challenging. This is true even when both 21-cm and optical 
velocities are known, for galaxy velocities are strongly clustered 
too. To be certain that every Parkes-detected source was correctly 
identified with its optical counterpart would require interferometric 
follow-up in every case - which is infeasible when one is dealing 
with hundreds of sources, as here. Wc have used interferometry, 
coupled with simulations, to be certain that only a handful of the 
remaining sources (less than 10) are still misidentified. This hand- 
ful is too small to affect the main correlations. 

The Equatorial Strip Survey (ES henceforth), reported on here, 
was a search through HIPASS cubes between declinations — 6°and 
-1-10°. Thus 5780 deg"^ of Hi data, in the velocity range between 
— 1280 and -1-12700 km s^^, were searched largely by eye to come 
up with 1107 sources (Garcia-Appadoo PhD Thesis 2005). The 
Equatorial Strip was chosen because: (i) it was approximately per- 
pendicular to the Galactic Plane and so it is mostly dark enough 
to make LSBGs detectable; (ii) it is accessible for follow up with 
a large range of instruments including the VLA, and for some of 
its area the SDSS-DR2; (iii) it includes the area searched for LS- 



BGs by Impey et al. (1996) using the APM machine (Cawson et 
al. 1987) to scan UK-Schmidt plates. This should eventually lead 
to an estimate of the number of LSBGs still missing from wide- 
scale optical surveys. So long as LSBGs remain a putative reser- 
voir for missing baryons (e.g. Fukugita ct al. 1998) it is important 
to know this. Of the ES area, 50 per cent will eventually be scanned 
by SDSS, and 35 per cent of it is already publicly released. We con- 
centrated our search in the 1700 deg^ of the ES released earlier in 
DR2 (Abazajian et al. 2004). 370 cross identifications were made, 
and then refined down to only 195 to reduce possible misidentifi- 
cations to a minimum. These 195 sources have SDSS-DR2 optical 
diameters (50 per cent light diameter and 90 per cent light diam- 
eter) between 0.16 and 6 arcmins (West et al. 2008). The SDSS 
pipeline photometry on sources of this size was known to be ex- 
tremely unreliable so we had to devise new techniques to reduce the 
data (West PhD Thesis 2005), which delayed the project by nearly 
two years. Eventually we will analyse many more ES sources al- 
ready observed with SDSS. But the results found for the first 195 
are clear and interesting enough to merit being published now. 

We have measured, at 21-cm, the peak flux, the integrated flux, 
the line width, and the spectral shape (e.g. two homed or Gaussian 
mainly), and in the optical the luminosities at u, g, r, i, z, two radii 
Rr)()(.g) and R9o(.9) containing respectively 50 and 90 per cent of 
the g-band light, the inclination, the morphology and a fibre nu- 
clear spectmm set either on the nucleus or the brightest HII re- 
gion. Given the strong correlations between some of the colours 
this amounts to about a dozen independent measurements, not far 
short of the 14 desired. Chiefly missing is a rotation curve to spell 
out the distribution of Dark Matter. Nevertheless we have sufficient 
information to hope that some important aspects of galaxy system- 
atics will emerge. With 13 properties there will be ~ (13x 12)/2 ~ 
80 possible correlations to look for. 

The previous largest survey of discs was reported in a remark- 
able, but largely vmremarked paper by Gavazzi et al. in 1996 enti- 
tled The phenomenology of disc galaxies. They fovmd: (a) the mass- 
to-light ratio in the NIR is virtually constant and equal to 4.6 in H- 
band solar units; (b) Population I indicators are all anti-correlated 
with mass; (c) Conversely the luminosity and surface brightness of 
old stellar populations all increase with mass; (d) Bulge compo- 
nents all increase non-linearly with mass; (d) All the above proper- 
ties are independent of either morphological type or environment. 

Previous blind Hi surveys, with follow-up optical or NIR data, 
have been carried out by Zwaan et al. (1997), Spitzak and Schnei- 
der (1998), Rosenberg and Schneider (2000), Minchin et al. (2003), 
Davies et al. (2004), Giovanelli et al. (2005). Some of the results 
are discussed in Rosenberg et al. (2005), Rosenberg and Schneider 
(2003) and Minchin et al. (2004). The main results can be sum- 
marised briefly as follows: 

(i) all but one galaxy (Minchin et al. 2005, 2007) seem to have 
optical counterparts; 

(ii) although spirals predominate, the counterparts cover a wide 
variety of morphological types and luminosities; 

(iii) there are significant, but not overwhelming, numbers of low 
surface brightness galaxies; 

(iv) to within the measurement errors all Hi galaxies have the 
same Hi column densities (~ W'^°-^^^°-^^ cm~^) and this is not a 
selection effect. 

The various surveys differ from one another partly in radio sensi- 
tivity, but mainly in the quality of their optical or NIR follow-up 
data. Our Equatorial Survey is distinct in two main respects: a) we 
have been fastidious in rejecting all sources where there might be 
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more than one galaxy in or near the beam; b) we have available high 
quality multi-band SDSS-DR2 data for every source, data which is 
both sensitive and uniform, across a wide dynamic range. Rosen- 
berg at al. (2005) have attempted a similar analysis using, instead 
of optical, 2MASS J-band data. 

Our most interesting results can be summarised as follows: 
a) All galaxies have the same mass-to-light ratio in H-band so- 
lar units. In other words Gavazzi's Law holds over a dynamic 
range of 8 mag in absolute luminosity even in Hl-selected galax- 
ies free of optical selection effects, b) There is a clear correla- 
tion between SB and Luminosity. Combined with (a) it reveals 
a correlation between dynamical mass and optical radius cubed, 
again over a large dynamic range, but now with more scatter. In 
other words there is roughly constant global-density for galaxies 
(~ 1 X 10"^'* gm cm~^). The Virial Theorem then implies the 
same angular velocity ~ 4 x 10~^^ radians s^^, implying they 
could have rotated no more than 40 times in a Hubble epoch, (c) 
All the galaxies have the same global surface density in Hi, that is 
to say the same Hl-mass divided by optical-radius squared (R^); (d) 
All the galaxies have exponential profiles throughout most of their 
extents so that Rgo is correlated very tightly with R^o ; (e) There is 
a clear colour-luminosity correlation in the sense that more lumi- 
nous galaxies are redder, which is well known in optically selected 
samples. 

The remainder of this paper is arranged by section as follows: 

(2) The Radio Obser-vations describe the Hi survey, the detec- 
tion process and the contents of the resulting catalogue. 

(3) The Optical Observations looks into the sensitivity of the 
SDSS for picking up large angular-size low surface brightness 
galaxies, and galaxies with low ( ^™ ) ratios. It then discusses 
the general problem of identifying reliable optical counterparts to 
Hi sources in blind Hi surveys. 

(4) The Diversity of Sources displays complete data sets for 
a handful of detected objects, ranging from giant spirals to in- 
choate low surface brightness galaxies, which represent important 
Hl-selected class types. A montage of the optical counterparts illus- 
trates their very wide variety. Finally the complete data set is laid 
out in the form of tables. 

(5) The Correlations in Galaxy Properties reveals the 5 sepa- 
rate correlations among their properties. 

(6) The Discussion argues that the correlations, despite the 
small size of the sample, are significant, not unduly affected by dis- 
tance uncertainties, and compatible with the TuUy-Fisher relation: 
Line- width ~ L", provided a ~ 1/3. If no more than 6 physical in- 
variants control galaxies, as we argue, then 5 correlations suggest a 
degree of organisation among gas-rich galaxies which is surprising, 
and which must strongly constrain theories as to their formation 
and evolution. 



2 THE RADIO OBSERVATIONS 

The observations were carried out with the Parkes 64-meter radio 
telescope using the Multibeam System (Staveley-Smith et al. 1996) 
in which 13 adjacent beams employing 26 receivers track the sky 
at a rate of 1 degree per minute, returning several times to the same 
piece of sky weeks apart to minimise interference. The main char- 
acteristics of the survey are listed in Table [T] 

The HIPASS catalogue, that is to say the pipeline catalogue of 
sources based on this raw data, has been released over the years us- 
ing techniques, and with results described in Barnes et al. (2001), 
Meyer et al. (2004), Zwaan et al. (2004) and Wong et al. (2005). 



Parameter 



HiPASS value 



Sky Coverage 
Integration time per beam 
Average FWHM 
Gridded FWHM 
Pixel size 
Velocity range 
Channel Separation 
Rms noise 
3(T Hi Mass Limit" 
Nh/ limit" 



<5 <+25° 

450 s 
14.3 arcmin 
15.5 arcmin 
4 arcmin 
-1280 to 12700 km s-i 
13.2kms-i 
13 mjy beam^^ 
lO'' d2 , Mb 
7.8 X lO^* cm-2 



" For AV= 100kms~i 
Table 1. Summary of HiPASS survey parameters 



Our Equatorial Survey, based on the same raw data, but covering 
only 5700 deg^ around the Equator (Garcia-Appadoo, PhD The- 
sis, 2005) was analysed much earlier (1999) with the intention of 
providing an early source list in time for SDSS-DR2 release. Al- 
though it used almost identical search methods to HIPASS, and 
indeed copies most of them directly, it appears marginally more 
sensitive than the HIPASS catalogue of the same region (Wong et 
al. 2005), perhaps because we had more time to follow up marginal 
detections with the more sensitive narrow-band system at Parkes 
(Zwaan et al. 2004) thus winning fainter signal from noise. The 
Equatorial Survey (ES hereafter) covers 14 per cent of the entire 
sky in a band right round the celestial equator from declination - 
6° to -1-10° . Of this area about 1700 deg^ is covered by the 2nd Data 
Release of the SDSS (Abazajain et al. 2004) which we will use as 
our source of optical information. Fig. [T] shows the distribution of 
sources around the sky, both those with SDSS-DR2 optical data and 
those without. The gaps caused by the Galactic Plane, as well as the 
strong peak at an R.A. of ~ 12 hours due to the southern extension 
of the Virgo cluster are easily visible, note that the source at 20h 
R.A. lies in an additional small area of SDSS-DR2 coverage that, 
due to its small size, is not shown in Fig.[T] Fig.|2]shows the distri- 
butions of ES sources in both peak and integrated flux. As can be 
seen the distributions for both the ES sample with (hashed) and the 
sample without (clear) optical data are very similar. This indicates 
that the optical sub-sample we will be analysing here is representa- 
tive of the sample as a whole and hence can be used to characterise 
the properties of Hi selected galaxies in general. 

The noise in the Parkes Multibeam data is complex, consist- 
ing of a mix of receiver noise, solar side-lobe emission, confu- 
sion and ripple in the 21 -cm spectra with a characteristic frequency 
of ~ 1200 km s~^ due to standing waves with the wavelength of 
52 meters set up between the dish surface and the prime focus 
cabin by continuum sources in the beam. With such complex noise 
the selection effects which enter into the discrimination of weak 
sources cannot be anticipated, and must be recognised retrospec- 
tively. Multibeam data is presented to the observer in cubes with 
axes in the R.A., declination and radial-velocity directions. The 
cubes are 8 degs on a side separated into 4 arcmin pixels and di- 
vided in the third dimension into 1024 velocity channels each 13.1 
km s~^ wide. Given that the original radio beam was 7 arcmin half- 
maximum in diameter the pixel signals are partially correlated, as 
are the velocity channels which have been Hanning-smoothed to 18 
km s~^. These cubes (102 of them in our case) were searched indi- 
vidually using a mixture of numerical algorithms and the eye-brain 
system. Algorithms are helpful, particularly in reducing labour; but 
none has so far proved anywhere as reliable or sensitive as the eye- 
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Right Ascension (hours) 

(a) Distribution of R.A.'s 




Declination [degrees) 

(b) Distribution of Declinations 

Figure 1. Spatial distribution of the ES sample, (a) shows the distribution 
in right ascension. The dashed lines indicate the range where there is no op- 
tical SDSS-DR2 data available, (b) shows the distribution of ES sources in 
declination. The unfilled histogram indicates those Hi sources without op- 
tical SDSS-DR2 data and the line-filled histogram the sources with optical 
SDSS-DR2 data. 



brain (e.g. Kilbom 2002). Each cube can be represented in 3 com- 
binations (q, 5), (a, V) and (5, V) where V is the radial velocity, 
and all 3 are searched before deciding a source is probably present. 
Such probable sources can be followed up by using the multibeam 
system in a more sensitive narrow-band mode. A combination of 
much higher velocity resolution (1.3 km s~^ per channel) and the 
ability to integrate indefinitely on a possible source, rather than for 
a total of 450 seconds in the scans, is a powerful filter against spuri- 
ous misidentifications (Zwaan et al. 2004). Note that the existence, 
or otherwise, of an optical counterpart is in no sense used in select- 
ing or rejecting a source, making the survey truly blind. 

The ES Hi sample was selected using methods almost iden- 
tical to HIPASS selection (Meyer et al. 2004, Zwaan et al. 2004). 
The noise is complex so the only way to establish the likely real- 
ity of sources is to follow up a sufficient proportion with the more 
sensitive 'Narrowband System'. The hundreds of putative HIPASS 
sources followed up gives confidence that significantly less than 
5% of the remaining ES sources could be spurious, far too small 
fraction to affect the correlations we are searching for [Section 5]. 

Fig. [3] shows the integrated flux as a function of the velocity 
width of the sources. It appears that detection by integrated flux 




og S 



(a) Distribution of peak fluxes (Sp^ak ) 




bg S. (Jy km s ') 

(b) Distribution of integrated fluxes (Si„t) 

Figure 2. Flux distributions for the ES sample. The unfilled histogram indi- 
cates those Hi sources without optical data and the line-filled histogram the 
sources with optical SDSS-DR2 data. 
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Velocity Width at 50% 



Figure 3. Selection limits in Velocity Width-integrated flux space. The the- 
oretical 3(7 limit for selection based on Sint (constant signal-to-noise) is 
shown by the dashed line and the 3cr limit for Speak selection ( av^q ^ 
shown by the solid line. Most en'or bars have been omitted for clarity. 
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Velocity Width (km/s) 



Figure 4. Selection limits in velocity width-peak flux space. The 3cr (39 
mJy) shown here (dot-dashed line) can be seen to be a good match to the 
selection limit of the data. The dashed vertical line shows AV50 = 3 chan- 
nels (39.6 km s~^). Most error bars have been omitted for clarity. 



is dependent on line width, with broad line sources being harder 
to find in the noise. In peak flux however it is much easier to set 
a clean selection criterion independent of line-width (see Fig. [JJl, 
as Kilbom (2002) and other searchers in the multibeam data have 
found. 

Fig. 5 shows flux limited curves N{S) ~ S-^/^ filled to peak 
fluxes (a) and integrated fluxes (b). Fig. |5(a)| is suggestive of a peak- 
flux completeness limit ~ 50 mJy. 

The raw velocities are corrected for solar motion relative to 
the Local Group. The resulting distribution of radial velocities is 
shown in Fig. (6] Approximate distances are derived in a similar 
way as in Koribalski et al. (2004), using D — vlg/Hq, where 
vlg ~ Vsys + 300 sin I cos b. Ho = 75 km s~^ Mpc~^ is as- 
sumed throughout (and where necessary a cosmology with Qm = 
0.3 and = 0.7). 

The Hi masses were derived from the integrated fluxes using 
fitting processes identical to HICAT and: 



Mhi{Mq) = 2.36 X 10^ X D(MpcY x / S{HI) dV 



(3) 



where D is the distance in Mpc and the integral J S{HI)dV is 
the integrated flux in Jy kms^^, equation (3) assumes that all 
the sources are optically thin and that the upper levels of the 21- 
cm transition are fully excited. Fig. |7] shows the distribution of 
Hi masses for the ES sample. 

Although in this paper we are not going to try and compen- 
sate for Hi selection effects we need to be aware of how they 
may have shaped our sample. The first worry concerns column- 
density-sensitivity. Will we have sufficient sensitivity to pick up 
LSBGs with correspondingly low Nhi (see equation (2))? It can 
be shown (Minchin et al. 2003) that the column-density-sensitivity, 
Nhi (cm~^) , of any survey can be worked out retrospectively from 
the fluxes F (Jy km s^^) and sizes 59 (Hi diameter in arcmin) of the 
faintest sources within it according to: 



atoms km 



where AV is their velocity width in kms^^. This is the 
Hi equivalent of equation (1) in the optical. For the ES it corre- 
sponds to ^ 10^^ cm~^, or a low SB limit, according to equation 
(2), where /I21 is the mean SB inside the 21-cm area of the galaxy: 




(a) Numbers of the ES sample by peak flux (Speafc) 




S.^^ {Jy km s-i) 

(b) Numbers of the ES sample by integrated flux (Sint) 

Figure 5. (a) Source count against peak flux. The histogram shows the num- 
bers found in each bin of peak flux, (b) the completeness of the ES sample 
in integrated flux. The histogram shows the numbers found in each bin of 
peak flux. The curve in both panels represents N(S) ~ 5^^/^ as expected 
for a flux-limited survey. 



M2i(B) = 27-2.51og 



iV^ 



1020. 



mag 



(5) 



For a typical (Mm / Lb) ~ 1 this corresponds to a central SB fio 
(B) ~ 25.7 or an effective SB /le// (B) ~ 27.5 which is very dim, 
dimmer than any optical catalogue covering a significant area (see 
Table[3]for definition of different SB's). Indeed it is so dim that one 
might question whether the SDSS could reach it. But here enters a 
surprising result recently discovered by Minchin using the Parkes 
multibeam to carry out a much deeper Hi survey (HIDEEP; 9000 s 
beam^^) of a small area, and with full optical follow up (Minchin 
et al. 2003). Despite its great sensitivity to low column densities 
(Nhi 2 x 10^** atoms cm~^) they found that there are no low col- 
umn density galaxies. Indeed they found that all Hi selected galax- 
ies have, to within the errors, the same column density [~ 10^"'^ 
cm^^ if spread over 5 effective radii, as in optically selected sam- 
ples (Salpeter & Hoffman 1996)], see Fig. [8] We find the identical 
uniformity of Nhi in the ES, see Fig. |l l(a)l 

Whatever the reason for Minchin's strange law (Minchin et al. 
2003), and it needs explaining, the dimmest LSBGs we can expect 
to encounter in the ES will have, according to equation (5), and the 
fact that Jl^ff =1^21 —2.1 mag for an exponential disc: 
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Table 2. Hi properties of the ES sample 



ES Name 


RA 


Decl. 


Speak 


Sint 


sys 


W20 


D 


log 




(J2000) 


(J2000) 


(Jy) 


(Jy km s~^) 


(kms-i) 


(kms-i) 


(Mpc) 


(Mq) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


HIPEQOO 14-00 


00 14 36 


-00 44 42 


0.075 


17.18 


3914 


290 


50.9 


10.02 


HIPEQ0027-01a 


00 27 47 


-01 09 39 


0.039 


6.60 


3848 


223 


54.2 


9.66 


HIPEQ0033-01 


00 33 22 


-01 07 01 


0.131 


17.24 


1972 


146 


30.1 


9.57 


HIPEQ0043-00 


00 43 31 


-00 06 49 


0.064 


13.81 


4124 


287 


60.6 


10.08 


HIPEQ0051-00 


00 51 57 


-00 28 25 


0.117 


14.84 


1616 


173 


26.0 


9.37 


HIPEQ0058-I-00 


00 58 50 


00 37 46 


0.048 


4.97 


5338 


156 


75.4 


9.82 



NOTE.- An extract of the table is shown here for guidance. It is presented in its entirety in the electronic edition of the Journal. 




2.0x10' 4.0x10' 6,0x10' S.Qxlo' 1.0 
Recessionol Velocity (km/s) 

(a) distribution of recessional velocities 




2-0 

log AV (km s"') 

(b) distribution of velocity widths 

Figure 6. (a) shows the recessional velocity distribution of the ES sample, 
(b) shows the distribution of the 50% velocity widths (AV50) for the ES 
sample. The unfilled histogram indicates those Hi sources without optical 
data and the line-filled histogram the sources with optical SDSS-DR2 data. 



fi,^ff{B){mag) = 23.9 + 2.5 log 



/ M{HI) \ 



(6) 



which should be comfortably accessible with SDSS for galaxies 
with {Mm / Lb)(3 < 5 (see Section 3). 

The Hi properties obtained for the ES sample are listed in Ta- 
ble[2] which is provided in full in the electronic edition of this jour- 
nal. The columns are as follows: 

Column (l).-ES source name. 



HI mass (Mj 

Figure 7. Distribution of Hi masses for the ES sample 




50 100 500 

AY (km s-1) 

Figure 8. The distribution of the Nj// column densities among the ~100 
HIDEEP galaxies observed by Minchin et al. (2003). The dashed line shows 
the sensitivity Umit. There appear to be no low column density objects 



Column (2) and (3).- Fitted Hi positions in right ascension 
and declination (J2000). 

Column (4).- Hi peak flux density, Speak- 

Column (5).- Integrated Hi flux density, Sint. 

Column (6).- Hi systematic velocity. Mays, measured at the 
20% level of peak flux density. 
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Figure 9. Distribution of velocity differences l\V/W20 between radio 
sources and their optical counterparts in the ES for the original 310 sources. 
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Figure 10. Distribution in positional differences between the radio and op- 
tical sources for the final list of 195 sources. 



Column (7).~ Velocity line width, W20, measured at 20% level 
of the peak flux. 

Column (8).- Distance, D, in Mpc (see text for description). 

Column (9).- The logarithm of the Hi mass, Mhi, calculated 
using equation (3) in units of solar masses. 



3 OPTICAL IDENTIFICATION AND OBSERVATIONS 

We determined that the best source of optical data to follow up 
the radio survey would be the SDSS. With its 5 photometric bands 
(u, g, r, i, z), it offers a unique and homogeneous data-base. The 
SDSS was however initially aimed at working on small faint ob- 
jects and the pipeline software was not designed to cope with the 
arcminute sized galaxies which turn up in the ES so that up to now 
most SDSS galaxy work has had to be confined to objects beyond 
z = 0.02 (e.g. Blanton et al. 2003). Identifying the various large 
galaxy problems, finding and validating solutions for them, and 
rewriting parts of the software has held up the ES for at least 2 
years. This is no place to discuss those corrections as an account 
appears in West et al. (2008). We summarise them only briefly 
to emphasise how necessary they are when using SDSS to study 
galaxies of arcminute size. 



3.1 Limiting surface brightness 

One legitimate concern about the SDSS, which is a relatively shal- 
low survey (~ 55 s exposure per point), is whether it will go deep 
enough to detect the kind of LSBGs one might hope to find in the 
ES. Tests show that the main source of noise in looking for dim 
extended objects in the SDSS is photon noise in the sky subtraction 
(West et al. 2008, Strateva et al. 2001). In that case it can be shown 
that: 



(7) 



where f(a) is the fraction of the sky level one can go down to 
and find objects of angular size 9 (diameter in arcsec) with a given 
signal-to-noise f(a) is that fraction for filter a, while N(a) is the to- 
tal number of photons accumulated from the sky arcsec"^ through 
filter a. For the (;-band SDSS N(g) ^ 120, and if one demands a 
S/N of at least 10 for detection: 



For mapped galaxies the 621 is typically 10 S50 where 6'5o is 
the half light radius (e.g. Salpeter & Hoffman 1996). As /i2i = Mo 
-I- 3.84, detection with SDSS in the g band requires: 



Ato(3) < 20.8 + 2.5logio (Oso) 



(9) 



In contrast LSBGs are generally defined to have j-io{B) < 23.0 
(see below). In other words, as you would expect of such a short 
exposure survey, the SDSS will be capable of picking up LSBGs 
only if they have large angular size (6^50 > 10 arcsec). In the case 
of the ES all the sources are nearby (^J 10,000 km s~^ away) and 
although some of our claimed identifications lie close to the SB de- 
tection limit (18), they appear visible because they have patches of 
light that are brighter than average, such as HIl regions. Of course 
we will not detect galaxies with relatively small fractions of Hi, 
dEs and dSphs for instance, unless the galaxies are very massive. 
We can quantify this limitation by combining (9) with (5) in which 
case, for typical values of log Nhi ^ 20.5 cm~^ (Minchin et al. 
2003) and sky brightness, /iskyig) = 22.5 g mag arcsec"^, the ES 
will only detect galaxies for which: 



{g)/{MH,/L 



> 15 



(10) 



7121(3) 5? 22.1 + 2.5 logxo 6 



(8) 



where 9 is the optical diameter (corresponding to 2 x 9go) in the g 
band. Gas poor galaxies will only be found if they are close by (i.e. 
9 large). In other words the ES may miss a significant fraction 
of light in the Universe coming from gas-poor galaxies. It is no 
coincidence that all the galaxies we do detect appear to be late- 
type. 

Many definitions of galaxy SB, and what constitutes low SB, 
appear in the literature. For instance there are two quite different 
definitions of effective SB that are current: which is the actual 
SB at the half-light radius, and /I50 which is the mean value of the 
SB within the half-light radius. To avoid even further confusion we 
list all these definitions together in Table|3]and relate them all to the 
central SB /io . The relationships are only valid for pure exponen- 
tial profiles. Where we have converted from SDSS colors to B-V 
system we have used Cross et al. (2004). 

What constitutes a low surface brightness galaxy must be 
a matter of convention. However most optically selected surveys 
show a distribution in SB, however defined, which is approxi- 
mately Gaussian with a FWHM ~ 3.0 magnitudes. Thus galaxies 
with a SB 1.5 or more magnitudes dimmer than the peak value 
could reasonably be defined as LSBG's and that seems to be a 
common, though not universal, convention (e.g. Impey & Bothun 
1997). For photographic, and for shallow CCD surveys, the peak in 
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Table 3. Definitions of surface brightness for galaxies with pure exponential profiles 



Surface Brightness measurement Definition 



MO 






Central Surface Brightness (SB) 


MSO (= 


MO 


+ 1.82™) 


The SB at the half-light radius (sometimes called 'effective SB') 


M50 (= 


Mo 


+ 1.12™) 


The mean SB within the half-light radius (also called 'effective SB') 


M90 (= 


MO 


+ 3.9™) 


The SB at the 90-per-cent-light radius 


M90 (= 


Mo 


+ 2.3™) 


The mean SB within the 90-per-cent-light radius 


V-21 (= 


MO 


+ 3.84™) 


The mean SB within the outermost 21-cm contour (see text) 



NOTE.- /JQ ill magnitudes corresponds to the central SB Iq where /(r) = Iqe and o is the scale length. 
Their total luminosity Lt = 2TTloa^\ R50 = 1.68 «, R90 = 2.32 R50 and (empirically) R21 ~ 5 R50 = 8.4 a. 



the disc SB distribution appears to lie close to the Freeman value 
of /io(B) = 21.65. Thus any disc dimmer than mo(B) = 23.0 fJ,(B) 
could reasonably be regarded as LSB. The median SB for the ES, 
allowing for average bulge-to-disk ratio of 0.1, is mo(S) ~ 22.5 
mag arcsec^, or about 0.8 mag arcsec^ dimmer than the Freeman 
value for optical surveys, and we find objects as dim as i-io{B) ~ 
24 mag arcsec^. 

3.2 Identification of optical counterparts 

The unambiguous identification of optical counterparts to 21-cm 
sources found in a blind survey like ours is by no means trivial. 
Because of the strong clustering of galaxies, both in angular and in 
redshift space, it is all too easy to find a plausible optical coun- 
terpart for virtually every 21-cm source. It is surprising, for in- 
stance, that in HOPCAT, the published optical catalogue for the 
4315 HIPASS sources in the Southern sky (Doyle et al. 2005) there 
is not one source, not one intergalactic gas cloud or dark galaxy 
without a plausible optical counterpart. We shall now estimate the 
rate of false identification, i.e. the probability of finding a random 
optical galaxy within a given distance, in both angular and redshift 
space, of any given Hi source. We shall assume, as the observations 
clearly suggest (Staveley-Smith et al. 1996) that optical galaxies 
and Hi sources are clustered together For an ES or HIPASS source 
the acceptable volume Wacc in which an optical counterpart could 
lie is a long thin cylinder, centred on the source, with its long axis, 
set by the radial-velocity uncertainties, along the line of sight. For 
an ES source at a typical radial velocity of 2000 km s^^, the an- 
gular uncertainties in position (= Ro) (up to 5 arcmin) correspond 
to ~ 50 kpc, while the velocity uncertainties AV, taking account 
of both radio and optical uncertainties, amount to Ho AV (~30 
km s~^) (= h) or half a Mpc. Given the correlation function: 

p{r)dV ^nodV[l + ^{r)] (II) 

where ^(r) — {-^) ^ * and no is the average number of plausible 
galaxies Mpc~^, it is possible to integrate the probability of finding 
a random galaxy within the volume "Vacc of the acceptable cylin- 
der To a very good approximation the number within a projected 
distance Ro (in Mpc) of the source is given by 

N{< Ro) ~ 1.8 no r^-" i?o ' (ro ~ 8Mpc) (12) 

Notice that the number is only weakly dependent on Ro 
(because of the strong correlation) and dependent on the radial- 
velocity uncertainty not at all. This last is counter-intuitive but 
arises from the long thin shape of the cylinder. The ends of the 
cylinder are so very far from the centre that finding highly corre- 
lated galaxies within the ends is very unlikely. Conversely, obtain- 
ing very accurate optical velocities for plausible galaxies in the field 



Table 4. Number of random galaxies within D (Mpc) A6' (rad) of radio 
position. 



Ro 


N« Ro) 


(1) 


(2) 


2 


78 


1 


35 


0.5 


13 


0.3 


8 


0.2 


5 


0.1 


2.1 


0.07 


1.3 


0.05 


0.8 


0.03 


0.4 


0.01 


0.1 



does not greatly enhance ones chances of making an unambigu- 
ous identification when the Hi velocity-uncertainties may still (for 
S/N reasons) be larger (and oddly enough, making blind Hi surveys 
with bigger dishes will not help because the characteristic sources 
will be proportionately farther away (Disney 2008). 

To turn equation (12) into numbers it is necessary to adopt an 
optical Luminosity Function for the putative galaxies. If we adopt: 

ip{M) = (/p.K"+'e-" (13) 

where k= iqOMm^-m)^ ^ ^ j ip, = 2.1 x 10"^ h~^ and 
M, = -20.04 - 5 log h (Blanton et al. 2003) then Table g] where 
Col. (1) gives the angular size distance in Mpc, Ro = -D AO, D 
is the source distance in Mpc, is the angular distance between 
the source and the optical counterpart in radians and Col. (2) is the 
number of random galaxies (down to 3 mag below M«) to be ex- 
pected within that distance. As an example, consider a source at a 
typical radial-velocity-distance of 2000 km s~^ in ES (~ 27 Mpc) 
where 1 arcmin corresponds to ~ 8 kpc. The typical positional un- 
certainty in HIPASS ~ 1.3 arcmin (Meyer et al. 2004, Zwaan et al. 
2004) corresponds to ~ O.OI Mpc. At that distance, according to 
Col. (2) of Table 4, there is a 10 per cent chance of finding a ran- 
dom galaxy of roughly the right radial velocity, i.e. clustered with 
the Hi source, within the positional uncertainty. And the chance of 
a misidentification could be higher still if one was prepared to con- 
sider, as plausible candidates, objects more than 3 mag below L, 
in the luminosity function (our assumption in table [3j- Note that 
Doyle et al. (2005) identify optical candidates up to 5 arcmin away, 
leading to the probability of a plausible misidentification at that 
distance of more than 1. 

To reduce these potentially serious misidentification problems 
in ES we threw out all (90 out of 310) sources which looked, on in- 
spection of the SDSS-DR2 fields, to have more than 1 plausible op- 
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Table 5. Petrosian Photometry of the ES sample 



ES Name 


RA 


Decl. 


u 


9 


r 


i 


z 


PetroRso 


PetroRgo 




(J2000) 


(J2000) 












(") 


(") 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


HIPEQ0014-00 


00 14 36 


-00 44 42 


13.77 ± 0.04 


13.13 ±0.02 


12.91 ±0.02 


12.82 ±0.02 


12.77 ±0.03 


21.4 ±0.4 


59.0 ±0.8 


HIPEQ0027-01a 


00 27 47 


-01 09 39 


14.72 ±0.04 


14.09 ±0.02 


13.85 ±0.02 


13.79 ±0.02 


13.87 ±0.04 


21.4 ±0.4 


43.6 ± 0.8 


HIPEQ0033-01 


00 33 22 


-01 07 01 


15.65 ±0.06 


14.86 ±0.02 


14.58 ±0.02 


14.44 ±0.02 


14.40 ±0.05 


20.2 ±0.4 


45.9 ± 1.6 


HIPEQ0()43-00 


00 43 31 


-00 06 49 


13.90 ± 0.03 


12.95 ± 0.02 


12.52 ±0.02 


12.31 ±0.02 


12.14 ±0.03 


13.9 ±0.4 


36.4 ±0.4 


HIPEQ0()51-00 


00 51 57 


-00 28 25 


15.42 ± 0.03 


14.60 ± 0.02 


14.23 ±0.02 


14.06 ±0.02 


13.99 ± 0.03 


6.3 ±0.4 


17.4 ±0.4 


HIPEQ0()58+()0 


00 58 50 


00 37 46 


14.90 ±0.03 


13.91 ±0.02 


13.43 ± 0.02 


13.19 ±0.02 


13.07 ±0.03 


10.7 ±0.4 


23.0 ±0.4 



NOTE.- An extract of the table is shown here for guidance. It is presented in its entirety in the electronic edition of the Journal. 



tical counterpart within the radio beam (FWHM ~ 14 arcmin) and 
we obtained accurate optical velocities for all the rest using either 
those provided in NED, the SDSS-DR2 fibre or, for 20 galaxies, the 
Dual Imaging Spectrograph on the 3.5 metre telescope at Apache 
Point Observatory. All candidates with optical velocities discrepant 
from the Hi value by more than half the 21 -cm line- width Al/20 
were discarded. In addition, a number of other sources were re- 
moved either because they extended across 2 or more SDSS-DR2 
fields, or because there was a saturated foreground star within 1 
arcmin. 

Altogether of the original 310 HIPASS sources with plausi- 
ble SDSS-DR2 galaxies within the HIPASS beam and at the right 
redshift (as defined above) 90 were thrown out because of multi- 
plicity i.e. for there being more than 1 (up to 5) good SDSS-DR2 
candidates, 20 were too extended and 5 too near a bright star. Fig. 
[9] shows the distribution of velocity differences AV/W20 for the 
original 310 candidates and Fig.[TO]shows the distribution in posi- 
tional differences between the radio and optical sources for the 195 
which remain in the final ES list. The tail in this plot is due to the 
Hi centroid not being centered in the optical source. 

Most of the identifications fall within 2 arcmin of the radio 
position, consistent with the errors in those positions measured in 
the general HIPASS catalogue using interferometry (Meyer et al. 
2004). Nevertheless there remain a tail of optical candidates up to 
7 arcmin (half the FWHM beam) away from the radio centroid. 
Because of clustering (see Table 3) we can not rule out the pos- 
sibility that a handful of sources (probably less than 10) remain 
misidentified. This is too small a number to invalidate the main re- 
sults. Nevertheless we should acknowledge biases in our sample 
against galaxies in tight groups (too many in the beam), galaxies 
that appear very large (overlapping SDSS-DR2 fields), and dark 
galaxies or intergalactic clouds for which we will too easily find 
plausible, optically bright alternatives. Most of these biases are dif- 
ficult to avoid and must exist to an equal or greater extent in other 
Hi selected blind samples. In particular the number of sources (~ 
30 per cent) discarded because of clustering within the beam high- 
lights the difficulty of measuring Hi mass functions which are not 
somehow adjusted for confusion. 

The data we used all came from SDSS Data Release 2 (Abaza- 
jian et al. 2004) using the pipeline from DR3 (Abazajian et al. 2005) 
but no SDSS catalogue data was used, for reasons outlined below. 
DR2 covers 3324 deg^ , about half of which overlaps the ES. SDSS 
pipeline photometry of large galaxies is very inaccurate for a num- 
ber of reasons. The pipeline shreds large galaxies into a number of 
pieces, circular apertures are not appropriate, the inclination is not 
properly taken account of and worst of all sky subtraction can sub- 
tract much of a large galaxy away from itself. Thus we had to find 



ways of measuring the significance of all these problems and de- 
vise alternative methods of handling the data. These are discussed 
at length in West (West PhD thesis) and are being published (West 
et al. 2008) so that the community wanting to use SDSS to work on 
nearby galaxies can make use of them. 

Table |5] lists the SDSS-DR2 optical data for the ES sample. 

Column ( 1 ).- ES source name. 

Columns (2) and (3).- Right ascension and declination 
(J2000) of the optical position of each SDSS-DR2 source. 

Columns (4) ?o fS).- Petrosian magnitudes in each SDSS band 
u,g,r,i,z, Unlike the native SDSS-DR2 Petrosian magnitudes, 
our Petrosian apertures are elliptical. The shape of each aperture 
is defined from two-dimensional Sersic fits and we use the g-band 
to measure both the shape and size of each Petrosian aperture. For 
some galaxies in the sample, a Sersic profile was not a good fit. In 
these cases, a circular aperture was used. 

Columns (9) and (lO).-The PetroR50 and PetroR90 values re- 
ported are the g-band semi-major axes that encompass 50 per cent 
and 90 per cent of the total Petrosian flux respectively. All of the 
SDSS data are described in West et al. (2008). 



4 THE DIVERSITY OF SOURCES 

In this section, we show just how diverse the galaxies in a blind 
Hi survey actually are. A quick way to grasp the wide diversity of 
types within the ES is to examine their SDSS-DR2 images, which 
can be seen online (http://www.astro.washington.edu/HIgalaxies) 
showing a g,r,i composite image. They range between red, early 
type barred spirals several times more luminous than L,, contain- 
ing proportionally small amounts of gas [(Mhi/Lb) < 0.2], and 
extremely low surface brightness objects hundreds of times fainter, 
very blue, consisting almost entirely of gas and so irregular in op- 
tical appearance as to be 'Inchoate', i.e. to be made of nothing 
more than irregular smudges with apparently unconnected higher 
surface brightness patches here and there. Given just how diffi- 
cult such Inchoates are to see at all on the SDSS-DR2 it is per- 
haps surprising that not a single source lacks an optical counterpart 
- however ragged it may be. There is even one elliptical galaxy, 
HIPEQ0I54— 00, which appears to have a thin disc and it is being 
investigated further. 

To consider the diversity of sources in more detail we con- 
sider four individual galaxies each representing a different char- 
acteristic 'class-type' of object that is quite common in the sam- 
ple. As a fiducial comparison one can pick say HIPEQ 1507-1-01 
(NGC5850) which is typical of the kind of galaxy that turned up 
in optically selected catalogues and which were afterwards exam- 
ined at 21-cm (e.g. Huchtmeier 1988). It is optically very luminous. 
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(c) HIPEQ0821-00 



(d) HIPEQl 145+02 



Figure 12. These images illustrate characteristic galaxy types within a blind Hi survey such as the ES, with their Hi spectra superposed, (a) shows a 'Hydrogen 
Giant'; (b) a massive LSBG; (c) an Irregular Galaxy and (d) an 'Inchoate Galaxy' - see text. 



1.7'" brighter than L,, has a low Mh//Ls =0.15 but nevertheless 
contains a significant amount of Hi (log Mh/ = 9.8) simply be- 
cause of its high luminosity. It is early type (T=3), red [{g — r) = 
0.63] and dynamically very massive (log Mdyn = 1 1.7 M0) where 
Mdyn = Rsoig) X iW2o/sin{i)f/G. 

The greatest number of galaxies in the ES are late type spirals 
a magnitude or so fainter than L, but gas-rich (i.e. median Mhi/Lb 
= 0.91 ±0.16) and blue. 

(i) Our first 'Hi Selected' type is a 'Hydrogen Giant', that is 
to say a galaxy with more than 10^° Mq of Hi. HIPEQ2036-04 
(NGC6941) contains 2 x 10^° M©, apparently an upper limit if 
one excludes other nearby companions - as here, and about 5 times 
the Hi mass of the Milky "Way (Binney & Merrifield 1998) [Fig 



|12(a)| . Their huge hydrogen content means that such galaxies can 
be detected far away and this one is at 6200 km s~^ . It is also a gi- 
ant optically, being ~ 3 L, in 3 but its gas-to-light ratio is MmfLs 
= 0.7 that is 5 times higher than the fiducial value for giants. Fig. 
|12(a)| shows an SABb (T=3) with a high surface brightness core 
but a low surface-brightness disc with widely spread spiral arms. 
The integrated colour is (g — r) = 0.79, while the fibre spectrum of 
the core is typical of an early type spiral. What is remarkable about 
these Hydrogen Giants is that they appear to have slowed their star- 
formation and so still contain as much mass in gas as in stars. The 
dynamical mass is high (log Mdyn = 11-3), but 1.5 times less than 
NGC5850, the fiducial galaxy. There are 20 in the survey, but given 
their generally high redshifts (median ~ 5500 km s~^) hydrogen 
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(a) Distribution of column densities 



10.5 



10 



S 



8.5 



* ! • 



1 10 
Optical radius (kpc) 
(b) The Hi mass/radius correlation 

Figure 11. (a) shows the Distribution of column densities for the ES sample 
where the densities are obtained by assuming 621 = 5 X 850 (opt) - see 
text, (b) shows the Hl-mass/radius correlation, i.e. constant column density. 
The dotted line shows the best linear fit to our data. 



giants cannot be common in space. The astrophysical challenge is 
to explain their delayed evolution from gas into stars. 

(ii) Our second class-type is the low surface brightness galaxy 
of the more massive kind - typically an anaemic spiral like 
HIPEQ1303+03 (UGC08153) shown in Fig. [T2(b)l (other good ex- 
amples are HIPEQ 1228-1-02 and 2337-1-00). It truly is low SB with 
a (S') ~ 24.3 and yet its dynamical mass is 10^° '-* Mq. For- 
mally speaking there are even more massive LSBGs in the sample 
but their discs are so extremely dim that estimating an accurate 
axis-ratio, necessary to calculate a dynamical mass, often becomes 
problematic. The total luminosity is 0.5 L, [Mg = -18.8] and the 
Hydrogen content is logio Mhi = 9.0 [Mh//Ls = 0.77]. The colour 
is blue, (g — r) = 0.27, as is commonly the case with LSBGs (Mc- 
Gaugh et al. 1994) but there are also some reddish ones in the sur- 
vey. The question of how significant LSBGs are in cosmic terms 
(e.g. Fukugita et al. 1998) hinges upon two difficult questions: on 
how massive they are and how common. HIPEQ 1303-1-03 and its 
like confirm, beyond question, that true LSBGs can be massive (see 
also Sprayberry et al. 1993). To what extent blind Hi surveys can 
compensate for the dramatic selection effects against LSBGs in op- 
tical catalogues is not yet clear. Certainly such surveys find healthy 



numbers of LSBGs. However Fig.ll3lshows that although N^r/ is 
more or less constant in the sample, there is a perceptible fall-off in 
the column density with SB, with a sharp cut-off at 10^*^'^ Hi atoms 
cm^^. This is due to the instrumental limit of the Multibeam sys- 
tem for 400 sec integration so it could still be missing some low 
SB objects (see Minchin et al. 2004). There are ~ 10 LSBGs with 
dynamical masses > 10^" M©, and as many again which are just 
as dim apart from a small bright core. 

(iii) Our third type is the irregular for which HIPEQ0821— 00 
(UGC04358) in Fig. |12(c)| serves as an example though they 
are very heterogeneous. ES irregulars are naturally gas-rich, 
HIPEQ0821-00 having (log Mm = 9.3) almost as much Hi as the 
Milky Way with an Mhi/Lb = 3.3. At Mg = - 16.4 it is ~ LJ20, 
has a very blue colour [(g — r) = 0.17] and moderately high SB 
CPsoid) ~ 22.5) though SB's among the type vary by 5 magnitudes 
and global SB's for such irregular objects are rather meaningless. 
There are at least 30 irregulars in the ES sample. 

(iv) The last 'type', the 'Inchoates', are so dim and faint that 
they could scarcely be found in any other but a blind Hi survey. 
Our example HIPEQ1145-F02 [Fig. 1 12(d)) can barely be seen on 
the SDSS, having a SB 7i5o(<7) = 24.8 at the very limit of the sur- 
vey. We detect it only because it has a lot of hydrogen for its lu- 
minosity (Mhi/Lb = 7) which at Mg = -14.31 is 1/200 L,. The 
'Inchoate' label for these objects derives from their apparent to- 
tal lack of organisation. More irregular than irregulars, they have 
no cores or obvious centres, and appear as merely haphazard en- 
hancements of SB at what appear to be HII regions. In addition to 
being extremely gas-rich (Mhi/Lb > 5) they are generally blue, 
thus HIPEQ 1 145-1-02, by no means extreme in colour, has B — V 
= 0.48 and (g — r) = 0.27. Other good examples of what is a vir- 
tually new type of galaxy are HIPEQ0238-I-00, 0240-1-01, 0958-1-01, 
1227-1-01 and 1256-1-03. We say 'virtually' because one of them, 
HIPEQ 1227-1-01, is the famous cloud serendipitously found by Gio- 
vanelli and Haynes (1989) and at first thought to be a protogalaxy. 
Indeed it would have been the easiest Inchoate to find as it has the 
highest Mhi/Lb (~ 22) of any object in the ES sample. 

Now that we have a dozen or so 'Inchoates' in the ES to study, 
some of the puzzling questions raised by the original Giovanelli 
and Haynes cloud return with even greater insistence. It is the com- 
bination of their properties which makes it difficult to explain In- 
choates (Salzer et al. 1991, Grossi et al. 2007). Their extraordinar- 
ily high gas-mass fraction indicates little integrated past star forma- 
tion, while their blue colours [as blue as (B— V) ^ 0.3] can only be 
explained with star formation that rises sharply to the present day 
(Bruzual and Chariot 2003). However, if such a rise is only a tem- 
porary burst, then the galaxies should soon fade by 1.5 magnitudes 
while reddening from (B— V) = 0.3 to 0.5, leaving behind an opti- 
cally undetectable dark Hi cloud (Leitherer et al. 1999). However 
there are no such dark clouds in our survey, and none among 4315 
HIPASS sources detected across half the sky (Doyle et al. 2005), 
all of which are optically detected in either SDSS or the SuperCos- 
mos Sky Survey (Hambly et al. 2001). Thus a bursting explanation 
for their blue colours seems less likely than either a truly young 
galaxy or a steadily rising star formation rate. However, GALEX 
measurements of 2 of them (1145-f02 and 1256-^03) imply current 
star formation rates that are only a factor of ~ 2 greater than the past 
average. And if they were truly young, why do not blind Hi surveys 
like HIPASS find their predecessor dark protoclouds? 'Inchoates' 
are indeed a puzzle. 

Whatever they are. Inchoates make a dramatic contrast with 
fiducial galaxies like NGC5850 and with the higher SB spirals 
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Figure 13. The surface brightness/column density coiTelation for the ES 
galaxies where SB=/I5q((;), see Table 3. 



which make up more than half the sources in the ES sample. To- 
gether with luminous gas-poor galaxies, hydrogen giants, massive 
LSBGs and irregulars they illustrate a diversity which suggests 
a wide scatter among several underlying fundamental parameters 
such as mass, age, angular momentum and binding energy. But 
such turns out not to be the case. 



5 CORRELATIONS AMONG GALAXY PROPERTIES 

We have measured 9 quantities in all, combined with the inclina- 
tion (b/a) which can be used to correct observed quantities to face 
on: mass of Hi, Hi line width (as defined in Section 2), luminosity 
(in g), four colours (it — g), (g — r), (r — i), and (i — z), and 
2 radii (R50 and Rgo in g band) all defined in Section 3. What 
we are chiefly missing is some kind of environmental parameter, 
which could be vital. However we have been so fastidious in ex- 
cluding closely packed galaxies that we cannot use environment 
for now. Later, with more interferometry of the dubious identifi- 
cations within groups, we can add that to the study. Even so it is 
worth pressing ahead without it, provided we remember the caveat, 
because Gavazzi et al. (1996) showed that for the late type galaxies 
which appear in the ES, environment appears to be unimportant. 

From the 9 basic measurements, and the inclination and as- 
sumed distance, one can construct other more familiar and per- 
haps more interesting parameters such as circular velocity (Vo), 
dynamical mass Mdyn Vq R50/G), surface brightness 7*50 (s) 
(~ Lg/R^o), column density (Nm ~ M{HI)/ R^q), brightness 
temperature Tm (~ Nm / ^W2o), angular momentum per unit 
mass (q ~ Vb-Rgo) and so on. With at least 15 such synthetic pa- 
rameters to work with there are ~ (15x 14/2) f» 100 correlations 
to look for We could have used some statistical technique but we 
elected to look at all the correlations by eye for now. Correlation 
coefficients, and their significance are given for each correlation. 
Thus one can discount obvious selection effects, and the already 
well-known correlations, to concentrate on looking for what might 
be new or surprising. 

We are of course far from the first to attempt such a system- 
atic search among galaxy properties. For instance Brosche (1973), 
Balkowski (1973), TuUy & Fisher (1981), Gavazzi et al. (1996), 
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Figure 14. (a) shows the relations between the {u — g) and (r — i) optical 
colours and (b) shows the relation between the {g — r) and {i — z) optical 
colours of the ES sample. 



Blanton et al. (2003) looked at optically selected samples but 
with added Hi data, whilst Minchin et al. (2003), Rosenberg et al. 
(2005), Kovac (2007) and Begum et al. (2006) have recently looked 
at Hi selected sets followed up either by optical or by NIR measure- 
ments from 2MASS. Elsewhere we shall compare results, but for 
now we examine only our own correlations with an eye unbiased 
by previous work. 

As a start it is important, given the rather small number in 
the sample (195), to establish that there is sufficient dynamic-range 
within its intrinsic properties to make correlation analyses worth- 
while. To look for relations between two galaxy properties y and x, 
of the form y = Ax"^, it is elementary to show that the error on m 
due to measurement errors Sy, Sx, as opposed to real scatter in the 
properties themselves, is given by: 



5m 
m 



1.5a 



jSy/y) 

A{magy) 



(Sx/x) 
A{magx) 



(14) 
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Table 6. Properties of the ES sample with optical SDSS-DR2 data. 



Property 


Median 


Min 


Max 


Range 


Relative error 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Mass Hi (Mq) 


109.38±0.40 


107.5 


IOIO.2 


6.8 


.03 


Luminosity (</) (Lq) 


J09.51±0.38 


107-2 


IOIO.6 


8.5 


.01 


Area (7r_R|p) (kpc^) 


10l-91±0.1 


10-0.2 


103.12 


8.3 


.01 


Line- width W20 (km s~^) 


102.12±1.78 


IOI.6O 


102.60 


2.5 


.09 


Surface brightness (Lqpc~^') 


J0l.85±0.07 


100-88 


102.99 


5.2 


.01 


Dynamical mass (Mq) 


J0l0.14±0.53 


108.5 


IOII.3 


7.0 


.06 


Colour (g — r) 


0.39 ±0.11 
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Figure 15. The (g — r) colour/luminosity relation. 

where magx = 2.5 logio{x) and the A{mag)'s give the rel- 
evant dynamic ranges. Table[6]sho-ws in Column (1) the most rele- 
vant properties measured for the galaxies in the ES sample. 

Column (2).- Median value of the measurement. 

Columns (3) and (4).~ The maxima and minima values respec- 
tively, where the few objects in the extreme wings of the distribu- 
tion have been ignored. 

Column (5).- Dynamic range in magnitudes. 

Column (6).- Shows the estimated relative measurement er- 
rors, 1.5{Sx/x)/ A{mag x). 

In most cases the relevant measurement errors are small enough (< 
0.1) not to dominate a correlation with another quantity in the table, 
the main exceptions being line-width W20 and colour. We do not 
use line-width except in the combination W20R/G, i.e. dynami- 
cal mass, where the relative measurement errors are reduced by the 
large dynamic range to < 0.06. 

First notice that the 4 colours are strongly correlated (Fig. 1 14b. 
as was to be expected, so there is considerably degeneracy among 
them and, as revealed already, to first order there is only 1 indepen- 
dent colour (Strateva et al. 2001). Many properties correlate with 
luminosity and Fig.[T5]shows the colour-luminosity diagram, which 
is well known for optically selected galaxies, i.e. more luminous 
galaxies are redder. When it comes to disentangling the physics go- 
ing on the actual slopes in all such correlation diagrams will be 
crucial. Thus Fig.[T5]implies that Lr ~ L^-^ roughly. The Pearson 
correlation coefficient is Vp ^ 0.61, corresponding to a probability 
P{r > Tp) > 99.9%. This is the first of the 5 correlations. 
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Figure 16. The surface brightness/luminosity relation in g band. The dotted 
line shows the best linear fit to the data. 



Fig. [16] exhibits the second correlation (r^ ~ —0.55, P{r > 
Tp) > 99.9%), i.e. between surface-brightness and luminosity and, 
not unexpectedly less luminous galaxies are dimmer. While this has 
been long suspected (e.g. Bingelli 1984) one could never be certain, 
in any optically selected sample, to what extent it was an artefact of 
the dramatic selection effects acting on surface-brightnesses which 
lie so close to the sky (e.g. Disney 1999). However here all the 
galaxies have been identified without regard to their optical proper- 
ties and it seems very likely that Fig.[T6]tells us something funda- 
mental about the relation between luminosity and radius with few 
selection effects. While the scatter is large the correlation is clear 
and very roughly speaking the surface brightness Ti(g) ~ L^g^ . 

The third correlation (rp ~ 0.75, P{r > Vp) > 99.9%), 
between Hi mass and optical radius, is illustrated in Fig. [TT| It is 
a tight one and suggest, that all galaxies have the same Hi column 
density, i.e. Mhi oc R^. This was discovered first by Haynes and 
Giovanelli (1984) in an optically selected sample, discovered again 
in an Hi selected sample by Minchin et al. (2003) and confirmed in 
an interferometer study by Rosenberg et al. (2005) in the Arecibo 
Dual Beam Survey and Kovac (2007) and our ES measurements 
concurs. For now constant surface density, is an intriguing puzzle 
which needs to be explained. 

The fourth correlation {Vp ~ 0.81, P(r > Vp) > 99.9%) we 
find is between luminosity {g) and dynamical mass shown in Fig. 
[TT] Over more than 3 orders of magnitude the g luminosity is tightly 
correlated with dynamical mass, a correlation that appears percep- 
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Figure 17. The dynamical mass/luminosity relation in g band. The dotted 
line shows the best linear fit to the data. 
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Figure 18. The Tully-Fisher relation in g band. 



tibly tighter than the Tully-Fisher correlation which we show in 
Fig.[T8](see later). The slope is close to 1, implying a direct propor- 
tion between dynamical mass and luminosity. This is not new either 
but it was discovered in the remarkable paper by Gavazzi et al. in 
1996. In order to compare our data with Gavazzi's, we transform 
ours into the iJ-band, which should be freer of dust and of tran- 
sient star-formation effects. 113 of our ES sources are in 2MASS 
and Fig. [19] illustrates the relation between their colour and their 
luminosity, which implies: 

logiLu) = (1.31 ± 0.05) X log{Lg) - (2.9 ± 0.4) (15) 

i.e. Lh ~ . If we use this to transform our Fig. 1 171 data into 
the H-band we reach Fig.|20|where our galaxies are plotted on top 
of the GOLDMine sample assembled by Gavazzi and colleagues 
(2003) for a mixture of bright discs and Virgo dwarfs. Over al- 
most 4 orders of magnitude, disc-like galaxies have Mdyn ~ Lh 
with a scatter of only ~ 1 mag about the relation in both axis 
directions. Whether they are selected optically (Gavazzi), or by 
their neutral hydrogen signal as here, the 635 galaxies all appear 
to have the same i/-band Mass-to-Light Ratio ~ 4 in solar units. 
The numerical value 4 is an artefact of using 7?5o((?) in the conven- 
tional definition of Mdyn, which we have so far adopted. However 
Mdyn oc R X {W2osin{i))^ and as the outermost HI test particles 
are typically at ~ 5 R^oig) in images 21 -cm galaxies (Salpeter & 
Hoffman 1996, Kovac 2007) a more accurate Mdyn would be ~ 5 
times greater than we and Gavazzi have implied. In that case the 
_ff-band mass-to-light ratio for both our sample, and the GOLD- 
Mine sample, is ~ 20, leaving plenty of scope for dark matter. 
Note that this ratio cannot be luminosity dependent (see Fig. 20). 
'Gavazzi's relation' as we call it, is a remarkably uniform one for 
any extragalactic data-set. 

The fifth and last correlation (r^ ~ 0.18, P{r > Vp) ~ 99.9%) 
we look at in Fig.[22]is that between the two separate radii measured 
by SDSS i.e, R50 and Rc,o- The very tight correlation, and the ratio 
between them of 2.3 means that despite the very wide range of 
galaxy types in the ES their outer light profiles are all well fitted by 
an exponential law and the ratio does not change with luminosity. 




Log L„o 

Figure 19. The (g — H) colour H band luminosity relation for the ES galax- 
ies detected by 2MASS. 



6 DISCUSSION 

We set out to look for correlations within a sample of galaxies free 
of optical selection effects. Apart from inclination and velocity we 
measured 9 parameters for each galaxy. Since the 4 colours appear 
degenerate we are left with 6 a priori independent observables, be- 
tween which we find the following 5 correlations, 4 old and 1 new 
(A) which are shown in TablelT] where the columns are as follows. 
Column (1).- Correlations observed. 

Columns (2) and (3).- Denote the distance dependance, d^, of 
the left and right-handed parameters in the correlation. 

Column (4).- Indicates the difference between the two expo- 
nents in Cols. (2) and (3) (see text for explanation). 

Column (5).- Shows the observables introduced into the cor- 
relation, that have not been used before. 

Column (5) in Table [7] is pertinent to the independence of 
the 5 correlations. Unless a priori reasons can be advanced for a 
dependance of the new observables on previous ones, then the 5 
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Figure 20. The dynamical mass-H band luminosity relation for our samples 
(circles) and for the GOLDMine sample (triangles). 



Table 7. CoiTelations in the ES sample. 





(1) 


(2) 


(3) 


(4) 


(5) 


A 


SB/Lum 


dO 


d? 


2 


SB, 690, d 


B 


Col/Lum 


dO 


d2 


2 


(9 - 1-) 


C 




dl 


d2 


1 


AV 


D 


R90/R50 


dl 


dl 





050 


E 


Mm/Rio 


d? 


d2 





S (21cm) 



correlations must be independent. It might, for instance, be argued 
that 6r,Q should depend on ^90, but that ignores the real differences 
in 'concentration' q = ^90/^50 between galaxies of different mor- 
phological type (e.g. q = 2.3 for pure exponential but q = 4.2 for 
de Vaucouleur profiles). Concentration is a crude but quantitative 
proxy for morphological type, and there is no a priori case for this 
to be the same among Hl-sources. Thus we can argue that the 5 
correlations are truly independent. 

Although the number of galaxies is modest the dynamic range 
of most of their properties is large (i.e. more than 7 magnitudes in 
luminosity, area, Mh/ and Mdyn) whilst the correlations are, in a 
statistical sense, highly significant. So far as the existence of the 
correlations (our main interest here) is concerned the sample-size 
is not an embarrassment, though larger samples will be needed to 
pin down the slopes of the correlations more precisely. There is 
moreover a great diversity of galaxy types in the sample, ranging 
between early-type giant spirals and tiny 'Inchoates' of such low 
surface brightness that they would usually be, and indeed were, 
missed in optical surveys. The difficulty of making a reliable census 
of galaxies in light alone is illustrated by Fig. [T6] where it can be 
seen that at all luminosities galaxies can vary in SB over 3 or 4 
magnitudes, making the dimmest of them, even the giants, hard to 
pick out from the background by optical means alone. The narrow 
SB range of earlier catalogues (Freeman 1970, Disney 1976) was 
almost certainly a selection artefact. 

Two correlations in particular need explaining. Correlation 
(E), i.e. M_H7 ~ could be related to the stability of a gas layer 
to gravitational collapse, and hence to a trigger for star formation 
(e.g. Kennicutt 1989, Schaye 2008). The other, the SB/Luminosity 




Log L„3 
(b) 

Figure 21. Top: The ratio M^y„/ R? ratio as a function of the g band lu- 
minosity. Bottom: The relation between the dynamical mass and R? for our 
sample (circles) and for the GOLDMine sample (triangles). 



correlation (A) has considerable scatter (Fig. II6I 1 making it diffi- 
cult to pin down an accurate value for /3 in E ~ , where {3 is 
not far from 1/2. Combining our data with the GOLDmine sam- 
ple, in order to increase sample-size and hopefully reduce scatter, 
yields E ~ L^^^ , ot T, ^ R, or Lh/R^ = constant or, consider- 
ing Gavazzi's Law, Mdy„/Rg = constant (~ 1 H atom cm"'' within 
the optical radius) which, see Fig. |21(b)] is a pretty good fit. So the 
mass - and luminosity - density of galaxies is more or less inde- 
pendent of their Luminosity. This is at least a better mnemonic for 
recalling the rather scattered SB/Lum law. If true the Virial Theo- 
rem implies that all HI galaxies spin at the same angular velocity 
and rotate ~ 40 times in a Hubble epoch. 

What about the Tully-Fisher correlation AV ~ L"? 
Gavazzi's Law (C): Lh ~ Mdyn = RgAV^ while above: Rg ~ 
Ly^ ^ , hence AT/ ~ . In other words our data are consistent 
with the TF law provided a = 1/3 [see also Fig. |18| . The problem 
with using TF directly is demonstrated in Table [6] where it can be 
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Figure 22. The concentration index of Hi galaxies as a function of lumi- 
nosity. A best fit line has a slope of only -0.017. The ratio should be 2.32 
for a perfect exponential disk. 



seen that the low dynamic range in AT/ will always make it diffi- 
cult to find the exponent a accurately, and it makes more sense to 
incorporate into the Dynamical Mass (= R AV^) instead, as 
it is in Gavazzi's Law (C) where it leads to both a clear correlation 
with Luminosity (Fig.l20t and perhaps a natural physical explana- 
tion ['all galaxies are made of the same stuff']. 

The distance to most galaxies in our radial velocity range can 
be problematic so it is important to ask if and how distance un- 
certainties will affect the correlations. Where there is a difference 
(last column) between the distance-dependant exponents x in our 
list one expects distance uncertainties to introduce extra scattering 
to correlations, as for instance in (A) but not in (E). It is interest- 
ing to note that the tightest correlations observed i.e. (D) and (E) 
are the ones unaffected by distance uncertainties, whilst the loosest, 
i.e. (A) [Fig. [T5l and (B) [Fig. 1161 should indeed be the strongest 
affected. This hints that distance uncertainties probably are increas- 
ing the scatter. On the other hand look at Gavazzi's relation (B) and 
in particular Figures[T7]and[20] There is no increase in scatter when 
we come to lower radial velocities where one might expect the rel- 
ative distance-uncertainties to have the most pernicious effect on 
the correlations. This suggests that although distance uncertainty 
can increase scatter, that increase must be rather limited here, and 
exculpates the uniform but rather crude distances we have adopted. 
Note also that only 1 observable is size, and hence distance depen- 
dant (R = dSgo) while 5 are not, i.e. Eopt, Nhi, {g — r), AV and 
Concentration Rgo/Rso- This implies strong correlations among 
intrinsic galaxy properties (e.g. colour and SB, AV and Nhi) i.e. 
ones that have nothing to do with either scale or distance. 

How constraining could the five discussed correlations (A) to 
(E) be on theories of galaxy formation and evolution? That depends 
on how many fundamental, independent invariant physical proper- 
ties galaxies might posses. At present it is hard to imagine more 
than 7: total mass, baryon fraction, age, angular momentum, ran- 
dom energy, radius and central condensation. Once a galaxy has 
virialised, the Virial Theorem will provide one relation between 
them, leaving only 6 independent properties. Ignoring interactions 
one might expect, to first order, that these properties might be con- 
served by each galaxy though, on a secular time-scale some ran- 
dom energy might be dissipated, while some baryons might be 



ejected. Thus five correlations (if they are truly independent) within 
a six-parameter set hint at a high degree of organisation among 
gas-rich galaxies, something partially foreshadowed by the early 
Hi pioneers [Brosche 1973, Balkowski 1973]. We are presently car- 
rying out a Principal Component Analysis to investigate the degree 
of organisation in this data. 

On the observational front work in progress by us includes 
enlarging the sample size as more SDSS and Hi data [e.g. from 
the Bonn 100 metre Multibeam Survey (Kerp et al. in preparation) 
becomes available. Our automated SDSS photometric pipeline for 
large galaxies should greatly accelerate progress (West et al. 2008). 
Hi interferometry is being obtained, first to sort out identifications 
in tight groups and second to measure more actual Hi radii to get 
a better understanding of the Minchin-Rosenberg effect. So far we 
have exploited neither the colours nor the fibre spectra but work is 
in hand to include them too (West et al. in preparation). Statistical 
work on mean cosmic properties, such as the mean Hi cosmic den- 
sity, depends on a clear understanding of ones selection effects and 
this is underway too (Garcia-Appadoo et al. in preparation). We are 
also working on the reverse List, that is to say a list of objects in 
the SDSS survey which we would have expected to contain mea- 
surable Hi, and thus show up in the Equatorial Survey, but which 
do not (West & Garcia-Appadoo in preparation). 
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